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Abstract 

We demonstrate that the concept of a conservation law can be naturally extended 
from deterministic to probabilistic cellular automata (PCA) rules. The local function 
for conservative PCA must satisfy conditions analogous to conservation conditions for 
deterministic cellular automata. Conservation condition for PCA can also be written 
in the form of a current conservation law. For deterministic nearest-neighbour CA 
the current can be computed exactly. Local structure approximation can partially 
predict the equilibrium current for non-deterministic cases. For linear segments of the 
fundamental diagram it actually produces exact results. 

1. Introduction 

Cellular automata (CA) are dynamical systems characterized by discreteness in space and 
time. In general, they can be viewed as cells in a regular lattice updated synchronously 
according to a local interaction rule, where the state of each cell is restricted to a finite set 
of allowed values. 

As in any other dynamical system, conservation laws play an important role in CA. 
Additive invariants in one-dimensional CA have been studied by Hattori and Takesue pQ. 
They obtained conditions which guarantee the existence of additive conserved quantities, 
and produced a table of additive invariants for Wolfram's elementary CA rules. The sim- 
plest of additive invariants, namely the number of active sites ("active" meaning non-zero), 
plays especially important role in CA dynamics. CA possessing such invariant, to be called 
"conservative CA", can be viewed as a system of interacting particles, as described in In 
a finite system, the flux or current of particles in equilibrium depends only on their density, 
which is invariant. The graph of the current as a function of density characterizes many fea- 
tures of the flow, and is therefore called the fundamental diagram. Fundamental diagrams 
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of conservative CA were investigated in [SJHj. For majority of conservative CA rules, fun- 
damental diagrams are piecewise-linear, usually possessing one or more "sharp corners" or 
singularities. Although the shapes of fundamental diagrams vary, there is a strong evidence 
of the universal behavior at singularities, as reported in jH] . 

Conservative CA appear in various applications, and some special cases have been studied 
extensively. Rule 184, which is a discrete version of the totally asymmetric exclusion process, 
is a prime example of such special case [3 13 El 13 UH EEB E2 EH1 El- Although much 
is known about this particular rule, dynamics of other conservative rules exhibits many 
features which are not fully understood, and more general results are just starting to appear. 
For example, M. Pivato ^Hj recently studied conservation laws in CA in a very general 
setting, deriving both theoretical consequences and practical tests for conservation laws, and 
provided a method for constructing all one- dimensional CA exhibiting a given conservation 
law. Another recent work 16 investigates universality and decidability of conservative CA. 
Unfortunately, there exists no general result explaining the shape of fundamental diagrams 
for conservative CA's, in spite of a remarkable progress reported recently for a special class 
of rules [T7] . 

In this paper, we will introduce a natural extension of the conservation condition to 
include probabilistic cellular automata (PCA). Deterministic conservative rules will then 
become a special case of conservative PCA. For the nearest-neighbour case, this allows to 
parameterize all conservative PCA by a set of three parameters. In the three-dimensional 
parameter space, nearest-neighbour conservative PCA are represented by a polyhedron, with 
deterministic rules located at its vertices. We will then show that the general shape of the 
fundamental diagram for nearest-neighbour PCA can be partially explained by a mean-field 
type approximation. 

2. Probabilistic boolean CA 

In what follows, we will assume that the dynamics takes place on a one- dimensional lattice 
of length L with periodic boundary conditions. Let Si(k) denote the state of the lattice 
site i at time k, where % G {0, 1, ... ,L — 1}, k G N. All operations on spatial indices % 
are assumed to be modulo L. We will further assume that Si(k) G {0, 1}, and we will say 
that the site i is occupied (empty) at time k if Si(k) = 1 (sj(fc) = 0). We will also define 
s(k) = {s (k), s^k), . . . , S£_i(fc)}. 

In a probabilistic cellular automaton, lattice sites simultaneously change states form 
to 1 or from 1 to with probabilities depending on states of local neighbours. A com- 
mon method for defining PCA is to specify a set of local transition probabilities. For 
example, in order to define a nearest-neighbour PCA one has to specify the probability 
w(si(k + l))|sj_i(A;), Si(k), s i+ x(k)) that the site Sj(fc) with nearest neighbors Sj_i(fc), s i+ i(k) 
changes its state to Si(k + 1) in a single time step. For the sake of illustration, consider a 
recently investigated PCA rule |18j . where empty sites become occupied with a probability 
proportional to the number of occupied sites in the neighborhood, and occupied sites become 
empty with a probability proportional to the number of empty sites in the neighborhood. 
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The following set of transition probabilities defines the aforementioned PCA rule: 



w(l|0,0,0) = w(l|0,0,l)=p 

w(l|0,l,0) = l-2p w(l|0,l,l) = 1 -p 

w(l|l,0,0)=p w(l|l,0,l) = 2p 

^(1|1,1,0) = 1 - P 10(1|1,1,1) = 1, (1) 

where p G [0,1/2]. The remaining eight transition probabilities can be obtained using 
w(0\a, b, c) = 1 — w(l\a, b, c) for a, b, c G {0, 1}. 

We will now proceed to construct a general definition of PCA, with arbitrary neighbour- 
hood size. Let n be a positive integer, and let C = {c\, C2, . . . , c n } be a set of integers such 
that Q+i = q + 1 for all i — 1 . . .n — 1. The set {si +Cl (k), Si +C2 (k), . . . , Si +Cn (k)} will be called 
a neighbourhood of the site Si(k). We will assume that the neighbourhood always includes 
the site Sj(fc), i.e., c\ < 0. 

Consider now a set of independent Boolean random variables X i v , where i — 0, 1, . . . L— 1 
and v G {0, l} n . Probability that the random variable X i)V takes the value o 6 {0,1} will 
be denoted by w(a\v), 

Pr(X iiV = a)= w(a\v). (2) 
Obviously, io(l|v) + w(0|v) = 1 for all v G {0, l} n . The update rule for PCA is defined by 

Si(k + 1) = -Xi,{ aifci (fc),« i+ca (fc),...,* i+0n (fc)}- (3) 

Note that {s(k) : k — 0, 1, 2, . . .} is a Markov stochastic process, and its states are binary 
sequences s(fc) G {0, 1} L . The probability of transition from x G {0, 1} L to y G {0, 1} L in 
one time step is given by 

L-l 

Pr(y|x) = Y[w(yi\{x i+C1 ,x i+C2 , . . . ,x i+Cn }). (4) 

i=0 

To illustrate the above formalism, let us consider again the rule (0) defined at the be- 
ginning of this section. For this rule, we have C = {—1,0,1}, and each lattice site i is 
associated with eight random variables A i>0 oo, -Xt,ooi, -Xi.oio, X ifi u, X itl(Xh ^,101, ^i.no, 
and Aj im . Probability distributions of these r.v. are determined by ([[J. For instance, 
w(l|l,0,0) = p, and therefore we have Pr(A ijl0 o = l) = p, and Pr(Aj 100 = 0) = 1 — p. 
Similarly, w(l|0, 0,0) = 0, hence Pr(A" ij000 = l) = and Pr(X i>000 = 0) = 1, meaning that 
Aj ooo is m this case a deterministic variable. 

3. Conservative rules 

If for any initial distribution \i 

^(E S ^)) =^|>(fc + l)) (5) 
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for all k = 0, 1, . . ., then the PCA defined by (|5J) will be called conservative. The expectation 
value of the sum Yl?iZo s i{k) can be interpreted as the expected number of sites occupied by 
"particles", therefore the condition (JHJ requires that the expected number of "particles" in 
the system is constant. 

Since the expected value of the random variable X i)V depends only on the vector v, we 
will now define function / : {0, l} n — > [0, 1] as 

f(x 1 ,x 2 , x n ) — E (X i>{xiiX2r __ >Xn} ) . (6) 

Proposition 1 The probabilistic CA rule defined by is conservative if and only if 

L-l L-l 

f( x i, x i+l, ■ ■ ■ , x i+n-l) = Xi (7) 
i=0 i=0 

for all {x , Xi, . . . , x L _i} G {0, 1} L . 

Proof. To prove that (JHJ implies (jJJ), it is enough to choose deterministic initial configura- 
tion, s (0) = x , si(0) = Xi, s L _i(0) = x L _i for any {x ,x 1 , . . . ,x L ^} G {0, 1} L . Then 

/L-l \ L-l 

and using (j3J) 

'L-l \ L-l 



. i=0 / i=0 



E » \ ^Sii 1 ) ) = Y1 E » ( X i,{^ +C1 (0)^ +C2 (0),...,s l+C , l (0)}) 
L-l L-l 

^ ] -fix (^.i^i + cj ,a:i + C2>---i :r i + c n }) ~~ ^ ] f( X i+Cl ] X i+C21 • • • J^i+Cn}) 

L-l 

= 5^/( a; i> a; i+i>--->a ; i+n-i)> ( 9 ) 



i=0 i=0 

L-l 



i=0 



where j = i + C\. Condition (jSJ) requires that £^ ^X^o* anc ^ fi« (Silo 1 mus t 

be equal, hence (J7J) follows. 

In order to prove that (J2J) implies we will use the definition of the expectation value 

(L-l \ L-l 

Y^s^k + l) = Yl J Pr(s(A; + l) = x)^x i , (10) 
i=o / xe{o,i} L i=0 



where x = {xq, x%, . . . , Probability Pr(s(k + 1) = x) can be written 



as 



Pr(s{k + 1) =x) = ^ Pr(x|y)Pr(s(£;) = y) (11) 
ye{o,i} i 



where Pr (x|y) denotes the probability of transition from y to x in a single time step. 
Combining the above two equations and changing the order of summation we obtain 



'2,-1 



,i=0 

L-l 

= X] Pr ( s ( fc ) = y) Pr ( x ly)$^ Xi 

xe{0,l} L ye{0,l} L i=0 

L-l 

= ]T Pr(s(fc)=y)]T ]T Pr(x| y )x,. (12) 

ye{0,l} L *=0 xG{0,l} L 

Note that the last sum, Xl x e{o i} L -P r ( x |y) x i> is simply equal to the expected value of Xi 
given the previous state of the system y. This expected value must be equal to 

E {Xi,{yi +C1 ,y i+C2 ,...,y i+Crl }) — f{Vi+cij Vi+cii • • • > Ui+Cn): (13) 

and as a consequence we have 

'L-l 



,i=0 

L-l 



£ Pr(s(A;)=y)^/( 2/ 

i+ci j Vi+C2 5 • • • 5 Vi+c„ ) 

y€{0,l} L »=0 

L-l 

= £ Pr(s(fc)=y)5> 

ye{o,i} L *=o 



precisely what we wanted to show. □ 
Following [HE], it is possible to obtain a simple condition which a probabilistic conser- 
vative CA rule must satisfy. 

Proposition 2 A probabilistic CA rule f is number- conserving if, and only if, for all 
{xi, X2, ■ ■ ■ , x n } G {0, l} n , it satisfies 

n— 1 

f(xx, x 2l ...,x n ) = x 1 + ^2 (/(°> 0, . . . , 0, x 2 , x 3 , . . . , x n ^ k+1 ) 

k=1 y k ' 

-/( 0,0, „ . ,0 ,xi,x 2 , . . . ,x n - k )). (15) 

k 

Proof. To prove that Condition (|T5jl is necessary, we will first note that /(0, 0, . . . , 0) = 0, 
which is a direct consequence of the condition (JJJ) applied to the configuration consisting of all 
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zeros, {xo, x\, . . . , xl-i} = {0, 0, ... , 0}. Consider now a configuration of length L > 2n — 1 
which is the concatenation of a sequence {x\,x 2 , . . . , x n } and a sequence of L — n zeros. If 
/ is conservative, such configuration must satisfy condition (J2J), hence 

/(0,0,...,0,x 1 ) + /(0,0,...,0,x 1 ,x 2 ) + --- 

+ f{x x ,x 2 , ...,x n ) + f(x 2 , x 3 , . . . , x n , 0) H 

+ f(x n ,0, . . . ,0) = xi + x 2 H \-x n , (16) 

where all the terms of the form f(0, 0, . . . , 0), which are equal to zero, have not been written. 
Replacing xi by in (fTBj) gives 

/(0, 0, . . . , 0, x 2 ) H h /(0, x 2 , . . . , x n ) 

+ f(x 2 , x 3 , . . . , x n , 0) H h /(x„, 0, . . . , 0) 

= x 2 H hv (17) 

Subtracting {HD from JEJ yields (fTKjl . 

To prove that condition (|T5|) is sufficient, we can apply it to each site x\ of a configuration 
{x ,xi, . . . ,x L _i}, 



f\%i-\-lj 2'j+2) ■ ■ ■ ) -^i+n) 

n— 1 

= X : 

k=l 



+i + (/(o,o, ...,o 

fc 

- /( 0,0,^. ■ ,0 , Xi+i, X i+2 , • • • , X i+n _ fc )) . (18) 

Obviously, when the above is summed over i e {0, 1, . . . , L — 1}, all the right-hand side terms 
except the first cancel, and we obtain (J7j). □ 
Local mappings / for conservative PCA exhibit some interesting symmetries. We have 
already demonstrated that /(0, 0, . . . , 0) = 0, and one can easily show that similar property 
is true when we replace zeros by ones, i.e., f (1,1, ... ,1) = 1. The following proposition 
shows yet another symmetry. 

Proposition 3 If a mapping f : {0, l} n — > [0, 1] represents conservative PCA rule, then 

J2 f(x 1 ,x 2 ,...,x n ) = 2 n ~\ (19) 

xi,a;2,...,a;ne{0,l} 

We will say that f satisfying the above condition is l-balanced. 

Proof. Consider a configuration t = {to, t\, . . . , ^l-i} G {0, 1} L , where L = 2 n . Now, let 
us construct a set of sequences of length n, A — {b^}^ ~q such that b^) = {tj, tj+i, . . . , t/ +n _i}. 
Superscript (j) denotes here just a consecutive number of the sequence b^ in the set A. 
Recall that all operations on subscript indices are taken modulo L. 
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Assume that we can find t such that all sequences in A are different. This means that 
each possible sequence of length n occurs in A once and only once, and therefore 

L-1 

^2 f(x u x 2 ,... 1 x n ) = ^2f(t i ,t i+u ...,t i+n _ 1 ) = N u (20) 

xi,a;2,...,a; n G{0,l} i=0 

where Ni is the number of l's in the sequence t. 

On the other hand, t G {0, 1} L , which is obtained from t by replacing all zeros by ones 
and vice versa, must also have the same property as t, i.e., 

L-1 

^2 f(xi,X 2 , ■■ .,X n ) = ^ f(ti,U+l,---,ti+n-l) =N , (21) 

asi,X2,...,x n e{0,l} *=0 

where N is the number of 0's in t. Comparing (|20|) and (|2*T|) we obtain N = N\ = 2"/2 = 
2™ _1 , exactly as required. 

The only problem left is to show that, indeed, for any n > 0, we can construct the 
sequence t of length 2 n (with periodic boundary conditions) such that all subsequences of 
length n occurring in t are different (and therefore, constitute a set of all possible sequences 
of length n). For example, for n = 3, t = 11101000 is the required sequence. One can see 
that sequences of length 3 occurring in t, 111, 110, 101, 010, 100, 000, 001, and 011, are all 
possible binary sequences of length 3. 

For a general n, the required t is equivalent to a hamiltonian cycle in the de Bruijn graph 
[IT?] of dimension n (or an eulerian cycle in the de Bruijn graph of dimension n — 1). It can 
be demonstrated that such a cycle always exists (more precisely, for a given n, there exist 
exactly 2 2 " ~ n of such cycles - see, for example, review article [20])- D 



4. Example: nearest neighbour conservative PCA 

We will now consider the case of nearest-neighbour PCA, for which the neighbourhood is 
defined by C = { — 1, 0, 1}. The most general PCA of this type is defined by eight parameters 
w(l\abc) G [0, 1] for a,b,c G {0, 1}: 

/(0, 0, 0) = w(l|000), /(l, 0, 0) = iu(l|100), 

/(0,0,1) = W (1|001), /(1,0,1) = W (1|101), 

/(0,1,0) = W (1|010), /(1,1,0) = W (1|110), 

/(0, 1, 1) = Ml|011), f(l, 1, 1) = w(l|lll). (22) 

Remaining eight probabilities can be determined by using consistency condition w(0\abc) + 
w(l\abc) = 1. Proposition |21 imposes extra conditions on /, 

f(x!,x 2 , x 3 ) = xi + /(0, x 2 , x 3 ) + /(0, 0, x 2 ) - /(0, xi, x 2 ) - /(0, 0, xi), (23) 

for all xi,X2,X3 G {0,1}. Writing these conditions explicitly, one obtains eight equations, 
among them only five are independent, implying that there are three free parameters in the 
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a 


P 


7 


Rule number 


f(x 1 ,x 2 ,x 3 ) 


1 








240 


Xi 


1 





-1 


184 


X\ — X\X2 + 











204 


X2 





1 


1 


226 


2 3 + XiX 2 - X 2 X 3 





1 





170 


X3 



Table 1: Five deterministic elementary conservative CA. 
solution. The general solution parameterized by a, (3 and 7 is given by 



/(0,0,0) 


= w 


;i|ooo) 


= 0, 


/(I, 0,0) 


= w 


;i|ioo) = 


a. 




/(o,o,i) 


= w 


;i|ooi) 




/(i, 0,1) 


= w 


;i|ioi) = 


a 4 


■P, 


/(0,1,0) 


= w 


;i|oio) 


= 1 — a — /3, 


/(i, 1,0) 


= w 


;i|no) = 


1 - 


P + l, 


/(o,i,i) 


= w 


;i|oii) 


= 1 - a -7, 


/(1, 1,1) 


= w 


;i|in) = 


1. 





where the three free parameters a, j3, 7 satisfy conditions 

ct>0, (3>0, a + (3<l, a + j<l, /3 - 7 < 1, (25) 

so that all probabilities w(l\abc) remain in the interval [0,1]. One can easily show that the 
following single expression for f(x%, X2, X3) is equivalent to (j23| : 

f{xx,x 2 , x 3 ) = 7(xix 2 - x 2 x 3 ) 4- axi + (I - a - f3)x 2 + (3x s , (26) 

where x%, x 2 , X3 G {0, 1}. 

The set of points (a, (3, 7) in 3D space satisfying (|25J) forms a polyhedron which has the 
shape of a pyramid with rhomboidal base and triangular sides, as shown in Figure There 
are five possible choices of parameters a, [3, 7 leading to purely deterministic rules (where 
all probabilities w(l\abc) are either zero or one), and they are listed in Table ^ These 
deterministic rules correspond to five vertices of the pyramid of Figure ^ 



5. Current conservation and fundamental diagrams 

In pp, Hattori and Takesue demonstrated that existence of a conserved quantity in a de- 
terministic CA is equivalent to a discrete version of a standard current conservation law 
dp/dt = —dj/dx. In what follows, we will show that a similar current conservation law 
holds for conservative PCA. 

We will first observe that the sum on the right hand side of the condition (J2J) can be split 
into two parts, one depending on x 2 , ■ ■ ■ , x n , and another which depends on xi, x 2 , . . . , x n -\. 
We will therefore define the function of n — 1 arguments 

n— 1 — ci 

J(x u x 2 , . . . ,x n _i) = - ^/(0,0, . . . ,0,xi,x 2 , • • -,x n _ k ) + y^ j x j , (27) 

*=i ' I ' j=i 

With this definition, it is straightforward to show that 
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Figure 1: Visual representation of the set of all conservative PCA rules in the parameter 
space (a, j3, 7). 



Proposition 4 A probabilistic CA rule f is number- conserving if, and only if, for all 
{xi, x 2 , ■ ■ ■ , x n } G {0, l} n , it satisfies 

f(x 1 ,x 2 , ...,x n )- x_ Cl+ i = J{xi,x 2 , . . . ,x n _i) - J(x 2 ,x 3 , . . .,x n ). (28) 

Applying this to all lattice sites we obtain 

f( s i+cii S i+C2i ■ ■ ■ 1 s i+c„) Si = ■/(Sj-f ci , Sj-|_ C2 , • • • , Sj_|_ Cn _ 1 ) 

— J( s i+c 2 i S i+c 3 ; ■ ■ ■ ) s i+c n -i)i (29) 

where we dropped k dependence for clarity. 

In order to understand the full meaning of the above equation, let us recall a general 
formulation of a conservation law in a continuous, one dimensional physical system. Let 
p(x,t) denote the density of some material at point x and time t, and let j(x,t) be the 
current (flux) of this material at point x and time t. A conservation law states that the rate 
of change of the total amount of material contained in a fixed domain is equal to the flux of 
that material across the surface of the domain. The differential form of this condition can 
be written as 

I -I- w 

Interpreting Si(k) as the density, the left hand side of (|29~j) is simply the expected change 
of density in a single time step, so that (j2Hj) is an obvious discrete analog of the current 
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conservation law (JHUj) with J playing the role of the current. We can pursue this analogy 
even further. For the partial differential equation (|30p one often assumes a functional relation 
between the current j and the density p, and we will do the same for conservative PCA: 
we will ask how the "average" current J depends on the "average" density in the stationary 
state, assuming that we start from a disordered initial condition. 

To be more precise, let us now assume that the initial distribution p is a Bernoulli 
distribution, i.e., at k — 0, all sites are independently occupied with probability p or 
empty with probability 1 — p, where p G [0,1]. Let us define p(i,k) = E^s^k)). Since 
the initial distribution is i- independent, we expect that p(i, k) also does not depend on i, 
and we will therefore define p(k) = p(i,k). Furthermore, for conservative PCA, p(k) is k- 
independent, so we define p = p{k). We will refer to p as the density of occupied sites. The 
expected value of the current J(s i+Cl (k), s i+C2 (k), . . . , Si+c^ik)) will also be i- independent, 
so we can define the expected current as 

j(k,p) = ( J (s i+ci (k) , s i+C2 (k),..., s i+Cn _ 1 (k))). (31) 

The graph of the equilibrium current j(p) = linn^oo j(k, p) versus the density p is known 
as the fundamental diagram. It has been numerically demonstrated [3] that for conservative 
deterministic CA the fundamental diagram usually develops a singularity as k — > oo, meaning 
that j(p) is not everywhere differentiable function of p. 

6. Current and fundamental diagrams for nearest-neighbour 
PCA 

For nearest-neighbour conservative PCA with the function / given by eq. (J2T?|) . the current 
becomes 

J (a?!, x 2 ) = Xx - /(0, Xt, x 2 ) - /(0, 0, x x ) = 'jx 1 x 2 + ax 1 - (3x 2 , (32) 
and the conservation condition (j2H| becomes 

f(xi,x 2 ,x 3 ) -x 2 = J(xi,x 2 ) - J(x 2 ,x 3 ). (33) 

For the five nearest-neighbour deterministic rules, the current J(xi, x 2 ) assumes particularly 
simple form. For rule 184, the expression (|3*2*|) reduces to J(xi, x 2 ) = xi(l — x 2 ), which means 
that the current J(si(k), Si+i(k)) is non-zero only if Si(k) = 1 and Si(k) = 0. This agrees 
with the interpretation of rule 184 as a system of interacting particles, in which a particle 
located at site i moves to site i + 1 only if the site i + 1 is empty. Rule 226 is similar, except 
that particles move in the opposite direction than in the rule 184. For rules 240, 170 and 
204, the current becomes, respectively, J(xi,x 2 ) = xi, J(xi,x 2 ) = —x 2l and J(x-i,x 2 ) = 0. 
This means that all particles move to the right (rule 170), to the left (rule 240), or stay in 
the same place (rule 204). 

It turns out that the simplicity of the expression for J(xi,x 2 ) for deterministic nearest- 
neighbour rules makes it possible to compute the expected current j(k,p) exactly. For rule 
204, we obviously have j(k,p) = 0. For rules 240 and 170 the expected current is just a 
linear function of p, and we have j(k, p) = p for rule 240 and j(k, p) = —p for rule 170. 
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Rule nr 
240 

184 

204 

226 

170 



J(x 1 ,x 2 ) 

Xi 

Xi(l - x 2 ) 


(1 - xi)x 2 
-x 2 



k+l 



_p ~^fttU+i 



i=l 



k+l 



j{k, p) 
P 

2k + 2 



k+l—il 



k+l+i 





2k + 2 



p k+1 -\l-p) 



k+l + i 



-P 



j(p) 



P '2 



1 
2 



Table 2: Expected current j(k,p) and the expected equilibrium current j(p) for five deter- 
ministic elementary conservative CA. 

For rule 184 (and its generalizations), one can also compute the expected current j(k,p) 
exactly, as done in [TT] . 



fc+i 



%=i 



2k + 2 



^P)-i-,-E^U + 7_ i j/ +1 -(i-/>) 

and one can prove that 



k+l+i 



(34) 



j(p) = lim j(k,p) 

fc^oo 



P ~2 



1 

+ 2- 



(35) 



This means that the graph of j(p) versus p has a singularity at p = 1/2, as shown in 
Figure Et- Rule 226 exhibits similar behavior, except that the direction of the current is 
reversed. Expressions for j(k,p) and j(p) for all five deterministic nearest-neighbour rules 
are summarized in Table 121 Unfortunately, these results cannot be easily generalized to 
compute expected currents for probabilistic rules. Even the equilibrium current j(p) for the 
general rule given by eq. (J26|) does not seem to be analytically tractable, not to mention 
the time-dependent current j(k,p). In the following section, we will use mean-field type 
technique to obtain approximate equilibrium current for probabilistic rules. 



7. Mean Field Current 

For nearest-neighbour rules, we can write the following expression for the average current, 
using equations (}3*Tj) and (J3*2~j) . 



j(k) = 7^ (U (s i _ 1 (A;)s i (A;)) + (a - 0)p, 



(36) 



where we used the fact that £ , /i (s i (/c)) = p. We will use abbreviated notation for probabilities 
Pr(sj_i(fc), Si{k) = axa 2 ), denoting them simply by Pfc(aia 2 ), where ai,a 2 G {0, 1}. -Pfc(ri) 
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is in this notation the probability that the pair (sj_i(A;), Si(k)) is at time step k in the state 
(1, 1), and since 

E^(s i ^ 1 (k)s i (k)) = Pr((s l ^ 1 (k),s l (k)) = (1,1)), (37) 

we obtain 

j(k)=-yP k (l,l) + (a-0)p. (38) 
By probabilities Pk{axQ>2) must satisfy the following relationship 

P k+ i(aia 2 ) = w{aia 2 \bib 2 b z h)P k {bib 2 hb / C), (39) 

&i,6 2 ,&3,&4e{o,i} 

where w(aia 2 \bib 2 b 3 b^) = w(ai\bib 2 b3)w(a 2 \b 2 b3b4), and the values of transition probabilities 
w are given by (J22j). As before, P k {bxb 2 b^b^) is the probability that the sequence of four 
consecutive sites at time k assumes value (bi, b 2 , 63, 64). 

The equation (JBTJj) gives two-site probabilities at time /c + 1 in terms of four-site probabil- 
ities at time k. In order to make it useful, we have to eliminate four-site probabilities from 
this equation. This can be done only approximately, using Bayesian extension formula also 
known as local structure theory 



P (h h h h \ P k(blb 2 )P k (b 2 b 3 )P k (b 3 h) 

Pk{bihhh) w — msm — • (40) 

where P k (l) = p, P k (0) = 1 - p. 

This approximation transforms (|39|) into a set of recurrence relations for pair probabilities 
Pfc(aid 2 ). Fixed point of this set corresponds to the stationary state, and the knowledge of 
the fixed point value of -Pfc(ll) allows us to find j(p) by ([3*8|). 

Although computation of the aforementioned fixed point is too long and tedious to be 
attempted by hand, it can be carried out by a computer algebra software. The resulting 
stationary value of the current is shown below (we used "mf" subscript to indicate mean-field 
nature of the approximation (|4"U1) ): 

j m{ (p) = (a-(3 + 1 )p- A 2t P (~ l \ 4- S ( 41 ) 

2( 7 p-7 + /5)(-7p + 7 + a) 



where 



A = (a - (3 + 7) (a - (3 + 7 - 7p) - a - (3 + a/3 



-sgn(7) a/co + cip + c 2 p 2 , (42) 

c = p 4 - 2 (2 7 + 1 + a) (3 3 + (-5a 2 + 4 7 - 4 7 a + 1 + 2 7 2 ) (3 2 
+2 (a - a 3 - 7 2 + 3«7 2 + 2 7 a 2 ) (3 - 2a^ 2 - 2a 3 

+4a 3 7 + 7 4 + a 2 - 4 7 a 2 + a 4 + 2 7 V, (43) 

d = -2 7 4 + 2 (a - 13) 7 3 + 2 (a + (3 - a 2 - (3 2 - 7a(3) 7 2 

+2 (/3 3 - a 3 + a 2 - (3 2 + 2/? 2 a - 2a 2 /5) 7 (44) 
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and 

c 2 = 7 2 (a 2 + f3 2 + 7 2 - 2 7 a + 6a/3 + 2 7 /?) . (45) 

The above approximation does a remarkably good job in predicting j(p) for deterministic 
rules. In fact, it yields all five expressions for the equilibrium current reported in the last 
column of Table 121 

For non-deterministic cases, the approximation does not work equally well. In order to 
evaluate its accuracy, we performed computer simulations in which the current j(k,p) was 
recorded for large values of k, thus approximating the equilibrium current j(p). Results are 
reported in Figure El where we show fundamental diagrams for a selection of five nearest- 
neighbour rules. Dotted lines represent simulation results, while continuous lines represent 
the approximate current given by (j41j) . Values of parameters a, (3 and 7 are shown in the 
upper left corner of each diagram. Only one diagram, namely Figure |2t, corresponds to a 
deterministic rule (rule 184). All other diagrams represent non-deterministic rules. 

Diagram (b) correspond to the rule located on the edge of the polyhedron of Figure Q 
joining rules 184 and 240, close to the vertex of rule 184. Similarly, diagram (c) represents 
the rule located on the edge joining rules 184 and 170, again close to the vertex of rule 
184. Remarkably, it appears that the local structure theory predicts correctly the shape of 
the linear part of the fundamental diagram, while the nonlinear part is only approximately 
correct. It also predicts existence of the singularity in the fundamental diagram, while 
numerical simulations indicate that no such singularity exists. 

When we consider a rule located on the edge joining rule 184 and rule 204, as shown in 
diagram (d), no part of the fundamental diagram is predicted correctly by the local structure 
theory (LST). The shape of the LST diagram, however, is not far from the observed shape. 
One can note similar phenomenon for a rule located on the base of the pyramid of Figure Q 
(diagram e) or in the interior of the pyramid (diagram f). 

The above evidence suggest some conjectures about fundamental diagrams of conservative 
PCA. First of all, singularity of the fundamental diagram appears to be a purely deterministic 
phenomenon. The fundamental diagram has a singularity (at p = 1/2) only for rules 184 and 
226. For all other rules, the graph of j(p) has no "sharp corners". Recall that the general 
probabilistic CA rule is defined in terms of probabilities w(a\v). If all of these probabilities 
belong to the set {0, 1}, such a rule will be called deterministic. All other PCA rules will be 
called truly probabilistic. 

Conjecture 1 The equilibrium current j(p) for any truly probabilistic conservative CA rule 
is a differentiable function of p for p e (0, 1). 

Secondly, linear parts of fundamental diagrams can be computed exactly by using local 
structure approximation. This has been observed in [Hj for deterministic rules, and appears 
to be true for probabilistic rules as well. 

Conjecture 2 If the equilibrium current j(p) of a probabilistic conservative CA rule is a 
linear function of p in some interval (a,b) C (0,1), then the local structure approximation 
for j(p) in this interval is exact. 
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Figure 2: Fundamental diagrams of selected nearest-neighbour conservative PCA rules. Hor- 
izontal axis represents the density p, while the vertical axis represents the stationary current 
j(p). Solid dots represent simulation results, and the continuous lines are graphs of Jmf{p) 
vs. p. 
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8. Conclusion 



We have demonstrated that the concept of a conservation law can be naturally extended 
from deterministic to probabilistic cellular automata rules. Local function / for conservative 
PCA must satisfy conditions analogous to conservation conditions for deterministic CA. 
We also proved that the local function must be 1-balanced. For nearest-neighbour rules, 
conservation conditions reduce the number of free parameters defining the rule to three, 
making it possible to visualize all nearest-neighbour rules as points of a polyhedron in 3D 
space, with deterministic rules located at vertices of the polyhedron. 

Conservation condition for PCA can be also written as a current conservation law. For 
deterministic nearest-neighbour rules the current can be computed exactly, but for proba- 
bilistic rules it appears to be a much more difficult problem. Nevertheless, local structure 
approximation can partially predict the equilibrium current. For linear segments of the 
fundamental diagram it actually produces exact results. 

Although our conjectures regarding fundamental diagrams are supported by numerical 
evidence for nearest-neighbour rules, obviously more evidence is needed, especially for larger 
neighbourhood sizes. Preliminary evidence indicates that both conjectures hold for four- 
input rules and five-input rules, and these findings will be reported elsewhere. 

Another interesting question is the issue of representation. In it has been demon- 
strated that conservative deterministic CA can be viewed as systems of interacting particles, 
or, in other words, that conservative CA rules can be defined using "motion representation" 
(as rules defining motion of individual particles). Similar concept has been introduced in 
|15j as a "displacement representation". As pointed out in [22], the motion or displacement 
representation is analogous to Lagrange description of the fluid flow in hydrodynamics, in 
which we observe a particle and follow its trajectory. Standard definition of the CA rule 
(as in eq. |3J), on the other hand, resembles Euler description of the flow, where the flow is 
observed at a fixed point in space and the dependent variable represents the amplitude of a 
field at that point. For deterministic conservative CA, both Langrange and Euler descrip- 
tions are possible. It seems that it should be possible to construct Lagrange description for 
probabilistic conservative CA as well, although it is not clear how to accomplish this task. 
The main problem is that the particles are not strictly conserved in PCA, only their expected 
number is. Since a given particle may disappear at any time, it is not possible to choose 
one particle and follow its trajectory indefinitely. Some averaging procedure will clearly be 
needed here. 
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